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Anisotropy effects on the finite-size critical behavior of a two-dimensional Ising model on a gen- 
eral triangular lattice in an infinite-strip geometry with periodic, antiperiodic, and free boundary 
conditions (be) in the finite direction are investigated. Exact results are obtained for the scaling 
functions of the finite-size contributions to the free energy density. With £> the largest and £< the 
smallest bulk correlation length at a given temperature near criticality, we find that the dependence 
of these functions on the ratio and on the angle parameterizing the orientation of the corre- 

lation volume is of geometric nature. Since the scaling functions are independent of the particular 
microscopic realization of the anisotropy within the two-dimensional Ising model, our results provide 
a limited verification of universality. We explain our observations by considering finite-size scaling 
of free energy densities of general weakly anisotropic models on a d-dimensional film, i.e., in an 
L x oo d_1 geometry, with be in the finite direction that are invariant under a shear transformation 
relating the anisotropic and isotropic cases. This allows us to relate free energy scaling functions 
in the presence of an anisotropy to those of the corresponding isotropic system. We interpret our 
results as a simple and transparent case of anisotropic universality, where, compared to the isotropic 
case, scaling functions depend additionally on the shape and orientation of the correlation volume. 
We conjecture that this universality extends to cases where the geometry and/or the be are not in- 
variant under the shear transformation and argue in favor of validity of two-scale factor universality 
for weakly anisotropic systems. 
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I. INTRODUCTION 

Bulk critical phenomena can be divided into distinct 
universality classes [1]. Critical exponents, certain crit- 
ical amplitude ratios, and the near-critical behavior of 
thermodynamic functions are identical for the members 
of such a universality class and are called universal, 
since they depend only on macroscopic properties such 
as near-critical correlation lengths and not on the mi- 
croscopic details of the system under consideration. For 
instance, the correlation length in the asymptotic critical 
domain, i.e., for asymptotically small positive or negative 
t = (T — T c )/T c , is, for isotropic systems, described by 

£ = £±,o|*r, T^T C , (1) 

where T c is the bulk critical temperature, v is a universal 
critical exponent, £±_o are nonuniversal critical ampli- 
tudes, and i?£ = £+,o/£-.o is a universal critical ampli- 
tude ratio. 

An extension to the concept of bulk universality con- 
cerns systems that are geometrically confined on a length 
scale L in one or more directions, where L is large com- 
pared to all microscopic length scales of the system, 
such as lattice spacings (see [2-4]; for reviews see, e.g., 
[5, 6]). Consider the free energy density / in units of 



k^T (this normalization is used without further mention- 
ing for all free energy densities throughout this work) of 
a d-dimensional system that is isotropic at such large 
distances. Assume that / may be split uniquely into a 
nonsingular and a singular contribution according to 

f(T,L)=f ns (T,L) + f s (t,L), (2) 

where the singular contribution / s is defined as that part 
of / that becomes singular in t at t = in the bulk 
limit L — > oo. If f s exhibits scaling, its behavior in the 
asymptotic critical domain of large L and small \t\, where 
Wegner corrections to scaling [7] are negligible, may be 
described by a scaling function T according to 

L d f s (t,L) =T(x), (3) 

with the scaling variable x = (L/£ !+ fi) 1 / v t. If T exists, 
it is expected to be universal. It describes the scaling 
behavior of the asymptotic singular part of the free en- 
ergy density for given system geometry and boundary 
conditions (be) for the bulk universality class under con- 
sideration for the isotropic case. 
With the bulk free energy density 

fb(T) = lim f(T,L), (4) 

L— >oo 

we may split the free energy density according to 
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f(T,L) = f h (T) + f ex (T,L), (5a) 
f cx (T, L) = f si (T, L) + f is (T, L), (5b) 
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where / ex is called the excess free energy density, 
f s [(T,L) = L _1 / s f(T) represents any surface or interface 
contributions, and /f s is called the finite-size contribution 
to the free energy density. Denote the singular parts of 
/ox, fsf, and /fe by / eXiS , / sf> , and / fSjS , respectively. If 
the leading contributions to these singular parts exhibit 
scaling, we may write in the asymptotic critical domain 

L d f CXtS (t,L) = A(x), (6a) 
L d f fStS (t,L)^G(x), (6b) 

with scaling functions A and Q. These, again, are ex- 
pected to be universal. The Casimir amplitude is defined 
as the critical value of G, i.e., by 

A = 5(0). (7) 

The physical motivation for investigating the scaling 
functions A and Q and the Casimir amplitude lies in 
their close connection to the critical Casimir force and 
its scaling function, see, e.g., the monographs [8, 9], and 
the recent overview article [10] (see also Appendix B). 
Rcf. [10] also provides a collection of recent theoretical 
and experimental results. 

Now consider weakly anisotropic systems, which are 
characterized by a single bulk correlation-length expo- 
nent v. We assume here that their bulk correlation 
lengths are related by a shear transformation to the 
isotropic case. Then v assumes in any direction the same 
value as for the isotropic case, which is supported by 
other investigations [11-15]. On the other hand, the £± ;0 
depend on the direction and are, for a (i-dimensional sys- 
tem, described by the surface of a d-dimcnsional ellipsoid. 
However, in any given direction, Eq. (1) holds and as- 
sumes its isotropic value. 

The situation is less clear for the scaling functions on 
the right hand sides of Eqs. (3) and (6), since they may 
additionally depend on the parameters describing the 
shape and the orientation of the correlation lengths ellip- 
soid, i.e., on d— 1 correlation length ratios and d{d—l)/2 
angles. These additional d(d+l) /2— 1 parameters may be 
organized into adxd symmetric matrix A with det A = 1 
[16]. The question arises whether such a function can be 
considered universal [17-19]. Dohm [18] interprets any 
dependence on A as nonuniversal. Diehl and Chamati 
[19] suggest to define universality only after transforming 
to an isotropic system by means of a shear transforma- 
tion. Here we advocate the interpretation that scaling 
functions even for weakly anisotropic systems are univer- 
sal, if they depend on the anisotropy only through its 
long-distance properties parameterized by A and not on 
any microscopic details of how it is realized. We will 
return to this issue in Sec. IV. 

For simplicity, we consider only systems, where the 
transformation to an isotropic system leaves the geom- 
etry and the be invariant. In particular, we compute 
free energy scaling functions for the specific case of an 
anisotropic two-dimensional Ising model on an infinite 
strip of width L and explain our explicit exact results by 



investigating systems confined to rf-dimcnsional films of 
width L. We consider periodic, antiperiodic, fixed, and 
free be in the direction of the width L of the film. For 
all cases we define a scaling variable by 

x^(L/^l)^t, (8) 

where £+q is the T > T c amplitude of the bulk correlation 
length £l perpendicular to the film boundaries with a 
corresponding asymptotic critical behavior 

a=d>r, T ^ T c (9) 

where tffi/Z-l = % 

We do not consider complications arising from strong 
anisotropics [20] , from subleading long-range interactions 
[19, 21], or from scaling violations for large x arising in a 
region of large L for fixed £l that manifest themselves in 
a nonuniform convergence of the leading singular part of 
free energy densities towards the respective scaling func- 
tion in the asymptotic critical domain [22] . 

This work is structured as follows. In Sec. II, we dis- 
cuss the two-dimensional Ising model on a triangular lat- 
tice in infinite-strip geometry. In Sec. II A, the model is 
set up and basic quantities are defined. In Sec. II B, we 
recall explicit results for the bulk correlation lengths and 
discuss their behavior near T c . In Sec. II C, we derive 
the scaling behavior of the singular contributions to / ex , 
f s {, and /f s for periodic, antiperiodic (Sec. II CI), and 
free (Sec. II C 2) be to the extent they are defined. We 
derive explicit expressions for the corresponding scaling 
functions and critical Casimir amplitudes. In Sec. Ill, 
we discuss the finite-size scaling behavior for free energy 
densities of d-dimensional anisotropic films with periodic 
and antiperiodic be (Sec. IIIB1) and fixed and free be 
(Sec. IIIB2) and compare our results to those found in 
Sec. IIC. The discussion of the results of Sees. II and III 
is delegated to Sec. IV. 

II. TWO-DIMENSIONAL ISING MODEL 

A. Basic definitions 

Consider a two-dimensional Ising model on a general 
triangular lattice with lattice constants di along the three 
lattice directions, see Fig. 1. For simplicity, we assume 
a ferromagnetic model with only nearest-neighbor cou- 
plings Ji along the sides of length di. Consider this model 
on an L\ x L 2 3 rectangular geometry, with the L\ direc- 
tion parallel to the "1" lattice direction. Let there be 
Ni layers in the "1" direction and N 2 3 layers in the "23" 
direction perpendicular to the "1" direction, so that 

Li = JVidi, (10a) 
L 23 = N 23 h u (10b) 

where h\ is the height of the elementary triangle with 
respect to the side of length d\, see Fig. 1. 
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L\ — L for per. and antiper. be 
L\ —> oo for ++, 00, and 0+ be 



FIG. 1: Triangular lattice with lattice constants di on an 
infinitely long strip of width L. 



Periodic or antiperiodic be are easiest implemented 
along one of the three lattice directions. We want this 
direction to be parallel to one of the rectangle edges and 
therefore choose the "1" direction. Periodic be are im- 
posed by identifying a line of spins along the "2" direction 
(or, equivalently, the "3" direction) with spins along such 
a line at a distance L\ in the "1" direction (one may think 
of the lattice as wrapped around a cylinder of circumfer- 
ence Li, whose axis points along the "23" direction). For 
antiperiodic be, in addition the signs of the couplings J\ 
along one side of one such line are reversed. 

Both free ("0") be, i.e., no further neighboring spin 
for the last spins on one edge of the rectangle, and fixed 
("+") be, i.e., a fictitious neighboring spin with fixed 
value +1 for each last spin on one edge of the rectangle, 
are easiest implemented if such an edge is along one of 
the lattice directions. For the rectangle defined above, 
this leads to the unique choice of imposing "00," "++," 
or "0+" be in the "23" direction, where the two entries 
refer to the two opposite sides of the rectangle. 

With /3 = l/(k B T) and K t = f}J u the Hamiltonian of 
this model reads 

m,n 

+ K 3 )• (ii) 

With the partition function 

Z(T,L 1 ,L 23 ) = Yl & ~ PH > ( 12 ) 
{* M =±i} 

the rectangle free energy density is given by 

/ rcct (T, Li, L 23 ) = -(LiL*,) -1 \nZ{T, L U L 23 ). (13) 



oo 




L 



FIG. 2: Bulk correlation lengths ellipse with largest and 
smallest correlation lengths £> and £< , respectively. The L di- 
rection is the direction perpendicular to the infinite direction 
of the strip and £l and £oo are the bulk correlation lengths in 
these directions. 

We are interested in the free energy density 

f f rec t(T, L, L||) for per., antiper. be, 
/(T, L) = lim < 

L ^°° [/ r ect(T,L||,L) for ++, 00, 0+ be, 

(14) 

of an infinitely long strip, i.e., for an L x oo geometry, 
with periodic, antiperiodic, 00, ++, or 0+ be in the L 
direction. For the width L of the strip holds 

{L\ for periodic and antiperiodic be, (15a) 
L 23 for ++, 00, and 0+ be, (15b) 
and the same relation between N, N%, and ./V23. 

B. Bulk correlation lengths 

For two-dimensional systems, the ellipsoid of Sec. I de- 
scribing the bulk correlation lengths in the asymptotic 
critical domain reduces to an ellipse, having a major ra- 
dius £> and a minor radius £< . Define a bulk correlation 
"volume" by 

Vco = (16) 

and an aspect ratio by 

r = £</£>, 0<r<l. (17) 

Let the direction of the major radius be rotated by an 
angle 6 with respect to the L direction of the infinitely 
long strip, see Fig. 2. 

Call £oo the bulk correlation length in the infinite- 
length direction of the strip and £l the bulk correlation 
length in the perpendicular direction. Basic geometric 
considerations provide the relations 

C 2 =^ 2 sin 2 + £< 2 cos 2 0, (18a) 

= ?> 2 cos2 6 + C< 2 sin2 0. (l 8b ) 
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between the various correlation lengths in the asymptotic 
critical domain. 

Let £j, i = 1, 2, 3 be the bulk correlation lengths in the 
ith lattice direction. It follows that 

f £l for periodic and antiperiodic be, (19a) 



for ++, 00, and 0+ be. 



(19b) 



Generalizing (1), the bulk correlation lengths behave, in 
the asymptotic critical domain, according to 



si — s±,orl 



(20) 



with £+o/£i^ = R^. For the two-dimensional Ising 
model, well-known exact results are 



i/ = l, 



i? e = 2. 



(21) 



According to Eq. (A22) of Ref. [23], the T > T c asymp- 
totic bulk correlation lengths are given by [24] 



where 



with 



7 (0 = 



= -ln 7 W, 

a« + VaW 2 - 46Wc« 
2cW : 



(22) 



(23) 



a« = 2^(1 + z|)(l + z 2 k ) + 4z jZk (l + zj), (24a) 
6W=^(1-^)(1-^), (24b) 

C W^(l_^ )(1 _^2 )5 (24c) 

with i, j, fc = 1,2,3 or cyclic permutations, and with 

Zj = tanhi'Q. (25) 

At the critical point the bulk correlation lengths diverge 
and therefore = 1. This condition may be written as 

z\ + z 2 + z 3 + z 2 z 3 + z 3 zi + ziz 2 - ziz 2 z 3 = 1, (26) 

so that 

1 — z 2 z 3 — z 2 — z 3 



z\. c {z 2l z 3 ) = 



1 - z 2 z 3 + z 2 + z 3 



(27) 



is the critical value of z\ for given z 2 and z 3 . Combin- 
ing (20)-(23), we may expand around the critical point 
through linear order in t resulting in 

z\ = z l c [z 2 ,Z 3 ) - — ■ ■ — "Try (28) 

(1 - z 2 z 3 + z 2 + z 3 ) 2 ^> Q 

and relate the & according to 



2(1- z 2 z 3 )(z 2 +z 3 ) {l-zl)(l+zl) (l-z|)(l+ 



(29) 



where terms of higher than linear order in t have been 
omitted. 

Geometric considerations, which we do not reproduce 
here, yield 



Ar(di,d 2 ,d 3 ) 
^1/6,^/6, d3/6) : 



(30) 



where the function 
1 
4 



(31) 



provides the area of a general triangle with side lengths 
t\, £ 2 , and l 3 . Combining Eqs. (29) and (30), and ob- 
serving that the area per lattice site A s i tc is twice the 
area of an elementary lattice triangle, i.e., 

^sitc = dih! = 2A tI (d 1 ,d 2 ,d 3 ), (32) 

we obtain, asymptotically close to T c , the relation 



V c , 



2(1 - z 2 z 3 ){z 2 + z 3 ) 



-) Ante, (33) 



(l-zf)(l-zf) \d, 
which will be needed below. 

C. Free energy 



In this section, we derive explicit results for the scaling 
behaviors of / ex , / s f, and / fs for periodic and antiperiodic 
be and of /f s for free be. We provide scaling functions 
for their contributions per bulk correlation volume, as 
well as standard scaling functions and critical Casimir 
amplitudes for them. The scaling functions and criti- 
cal amplitudes will be expressed in terms of the scaling 
variable x from (8) and the anisotropy parameters r and 
8. For further reference we note that with (8) and (21) 
follows 



for T > T c , 



Lt/(M = x 
-Li/d L o = -2i for T < T c . 



(34) 



For the two-dimensional Ising model, the leading sin- 
gular behavior of / in the bulk limit L — > oo behaves 
as oc t 2 ln|£|. The arbitrariness in splitting the constant 
under the logarithm between singular and nonsingular 
contributions to / prevents the required unique splitting 
of / according to (2) and causes a violation of scaling. 
Consequently, T does not exist. On the other hand, such 
a splitting of /f s and therefore the function Q exist for all 
cases considered in this work. The existence of a corre- 
sponding splitting for / ex and / s f, and thus the existence 
of A, depends on the be, as we will see below. 

It is useful to define the strip free energies per site / 
and per bulk correlation volume / by 



f(T,L)=A site f(T,L), 
f(T,L) = V co f(T,L), 



(35a) 
(35b) 
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respectively. Let analogous definitions hold for /b, / ox , 
/ s f, /f s , and, if a unique separation as in (2) is defined, 
also for their singular parts. 



1. Periodic and antiperiodic be 

The free energy per site defined in (35a) is, according 
to Eq. (71) of Ref. [25], given by 



f(T,L) = -C + f sf (T,L) 

N-l 



V l /i(^) + V/i 2 (^)-/ : 



with 



2N 

3=0 

s 2^_+ I) 
3 N 
, _ 2nj 
3 = ~N~' 



(36) 



for periodic be, (37a) 
for antiperiodic be, (37b) 



and where [26] 

h{4>) = A -A lC oscp, (38a) 
/ 2 (0) = yj (A 2 + As) 2 - 4A 2 A 3 sin 2 (0/2), (38b) 



A 
A 3 



(1 + z\){l + z 2 )(l + z 2 ) + 8z lZ2 z 3 , 
2z 1 (l-z 2 )(l-z 3 2 ), 
2z 2 (f-z 2 )(l-z 2 ), 
2z 3 (f-z 2 )(l-z 2 ), 



and 



C = In 



(39a) 
(39b) 
(39c) 
(39d) 



(40) 



with the qualification that the interface function / s f is 
missing in Ref. [25]. For T > T c indeed / s f = 0. For 
finite L, there is no phase transition and therefore the 
free energy per site for T < T c is obtained by analytic 
continuation of (36). While for periodic be, this leaves 
/ s f at zero, the j = term has to be treated separately 
for antiperiodic be and we obtain 

/ s ( f p) (T,L)=0, 



(41a) 

f^(T, L) = ®£) ln (l-Zi)(l-^ 3 ) (41b) 



N 



(l + z 1 )(z 2 +z 3 ) ' 



with the Heaviside step function 9. The interface contri- 
bution for antiperiodic be for T < T c was overlooked in 



Ref. [25], leading to erroneous results there for T < T c . 
In particular, the 9 < plots for R = 00 in Fig. 3(a) and 
for L y = 00 in Fig. 3(b) are missing a linearly rising part 
towards 9 — > —00. This leads to incorrect statements 
between Eqs. (78) and (79) of Ref. [25] about the 9 < 
behavior of the infinitely long cylinder. 

Taking the limit L — > 00 of (36), we obtain the bulk 
free energy per site as 



-c 



^ Jo 2 



(42) 



and we note that the critical point corresponds to 

A - At - A 2 - A 3 = 0, (43) 

which is equivalent to (26). 

Combining (36) and (42) according to (5), we obtain 
in an obvious notation 



k(T,L) 



In 



From (38) we obtain 



1 ,2tt 1 n-i 



(44) 



f?(<P)-fl(<f>) = 

(A - Ai + A 2 + A 3 ){A Q -A 1 -A 2 - A 3 ) 

+ 4 [A A 1 + A 2 A 3 - A\ cos 2 (0/2)] sin 2 (0/2). (45) 

Near T c the quantity A — A\ — A 2 — A 3 is small and an 
appropriate approximation for < < 7r may be written 

as 



ff(<f>) - /! (0) ~ 2 A / A A, + A 2 A 3 - A\ cos2 J sin J 



^(Ao-AM^A^ ^-^_ n ^ (4g) 
where we have used that through linear order in t 



' (A - Ai + A 2 + i4 3 )(i4 - A x - A 2 - A 3 ) 



{Ao- A 1 )A 1 + A 2 A 3 



Lt/^ =x for T^T C 



(47) 



which follows from (27), (28), (34), and (39). An approx- 
imation analogous to (46) holds for tt < < 2ir, which is 
found by replacing — > 27r — there, under which (45) 
is invariant. Thus we may write 
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MT,L) 



1 

4tt Jo 



1 

2N 



N-l 

E 



In 



A - A x cos <f> {j) + 2^A a Ax + A 2 A 3 - A\ cos 2 %i sin 



+ 



(A - A X )N 



1 /-27T -, JV-1 

JU j=0 



(48) 



Using the results (Ala) and (A2a) for periodic be and 
(Alb) and (A2b) for antiperiodic be, we obtain 

N*f^\T,L) « 

(49) 



with 



(z) = -I d W In (l ± e -^+^) , 



(50) 



and where we note that 

I<+>(0) = -7r/12, =tt/6. (51) 

Close to the critical point, we obtain, after some alge- 
bra and with the help of (33), 



^/(A -A 1 )A 1 + A 2 A 3 {Zi\ 2 A 



(A - A\) ' \d x J V co 

Combining this with (49) gives, close to T c 



(52) 



^ /a) (T,L)«||^/(+/-)(5). (53) 

Using (19a), (28), (33), and (34), we obtain for the in- 
terface contribution (41b) to the free energy per site for 
antiperiodic be for temperatures close to T c 



&(T,L)^-^ 2 ^e(-x)x. 



(54) 



For periodic and antiperiodic be, the free energy con- 
tributions / cx , f s f, and /f s are dominated by their lead- 
ing singular parts / eX;S , / sf;S , and / feiS , respectively. With 
(19a), (34), and (35), we obtain the asymptotic singular 
finite-size and interface parts of the free energy per bulk 
correlation volume as 



'/(+/-) (x)/x 2 for T>T C , 
&\t,L)={ (55) 
'/(+/-) (x)/{Ax 2 ) for T<T C , 



and 



f^(t,L) = -Q(-x)/(4x), 



(56a) 
(56b) 



respectively. With (5b) follows the asymptotic singular 
excess free energy per correlation volume 



f!&{t,L) = fi»}(t,L), 



for T > T c , 



[I { ~\x) - x]/(4x 2 ) for T < T c . 



(57a) 



(57b) 



Eqs. (55)-(57) represent central results of the current sec- 
tion. Their most remarkable feature is that their right 
hand sides depend only on the scaling variable x and 
not on the parameters r and 9 describing the anisotropy 
of the system at large distances (nor on any other de- 
tails of how the anisotropy is realized). We conjecture 
that the right hand sides of (55)-(57) represent univer- 
sal scaling functions in the bulk universality class of the 
two-dimensional Ising model. In Sec. IIIB 1, we will ex- 
plain their independence of r and 9 in the more general 
context of d-dimensional films. 

Next we derive standard free energy scaling functions 
and critical amplitudes. Combining (16) and (18b) gives 



Vco/& = r cos 2 9 + r' 1 sm 2 



(58) 



so that multiplying (55) and (56) by L 2 /V co , while ob- 
serving (34), gives asymptotically 



L 2 f^ a) (t,L) = g^\x,r,6) 
/(+/") {£) 



r cos 2 9 + r 1 sin 2 9 ' 



(59) 



and 



L*f£(t,L) = - 



(60) 



Q{-x)x 
r cos 2 6 + r _1 sin 2 0' 

respectively. Therefore, we have 

L 2 f^(t,L)=A^(x,r,0) 

'Q^ (x, r, 9) for periodic be, 

l(-\x)-Q{-x)x £ (61) 
for antiperiodic be, 



■ r cos 2 9 + r 1 sin 2 8 



where _4,( p / a ) (x, r, 9) and Q^/^ [x, r, 9) are the scaling 
functions of the excess free energy density and the finite- 
size contribution to the free energy density, respectively. 
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0.4- 



which wc expect to be universal functions of r and 9 
within the bulk universality class of the two-dimensional 
Ising model. 




FIG. 3: Scaling functions g£>(x) (solid, < 0), G^(x) (solid, 
> 0), 0i ( ° o O) Or) (dashed, < 0), g£ +) (x) (dot-dashed, < 0), 



and gZ '(x) (dashed, > 0). 



For the two-dimensional Ising model, they generalize the 
isotopic-case structure (6) to the anisotropic case. We 
conjecture that _4(p/ a ) and Cj(p/ a ) are universal functions 
of x, r, and 9 for all infinite-strip systems in the bulk uni- 
versality class of the two-dimensional Ising model with 
the appropriate be in the L direction. Note that the 
scaling functions depend on the anisotropy parameters r 
and 9 only through a geometric factor. 

For the isotropic case r = 1, the results (59) and (61) 
reduce to 



and 



iL p o /a) (*) = 



^ o /a) 02) 



/(+/") (x) (62) 
for periodic be, 



\x) — 0(— x)x for antiperiodic be, 

(63) 

and thus we may relate 

G (p/a) (5;,r,9) = (r cos 2 9+r' 1 sin 2 6>) _1 G ( Z /a) (x), (64a) 
A ip/a) (£,r,9) = (r cos 2 9+r- 1 sin 2 9y 1 A^ a) (£). (64b) 

With x± — x, our result for Gfa£(x) is identical to 
6_ L (a;_ L ,0) from Eq. (57) of Ref. [27]. Even for the 
isotropic limit, and A^ appear not to have been 
previously published. (x) and G^o 0*0 are shown in 
Fig. 3. A cross check of G^J^ and A^^ with published 
scaling functions for the Casimir force is provided in Ap- 
pendix B 1. 

With (51), the critical Casimir amplitudes A(r, 6) = 
G(0,r,9) corresponding to (59) are 



A< P >M) 
A^(r,9) 



-7i712 



r cos 2 9 + r 1 sin 2 9 ' 
tt/6 

r cos 2 9 + r -1 sin 2 9 ' 



(65a) 
(65b) 



2. Free be 

The free energy per site defined in (35a) is, according 
to Eq. (A15) of Ref. [23], for free be given by [28] 



(7, L h L) = - ln(2 cosh K x cosh K 2 cosh K 3 ) 

^lll|J.-^ 1 c" ? | 



N 



+ ^ln(cosh J fC 2 coshi^ 3 ) ~^J2 ln ^' 

II # 

^^In^A^+p-A^ 1 ), (66) 



with 



•d = 



(2p - 1)tt 

iVii ' 



p= 1,...,N« 



P± = ± 
A± 



A-\ T ~a\E\ 



V /(A-F)2+4|SC|' 
A + F± V /(A-F) 2 + 4|SC| 



and 



The strip free energy per site is obtained as 
P°\T,L)= lim /^(T,L||,L) 

L||— s-oo 

= — In (2 cosh i^i cosh K 2 cosh i^) 



+ — ln(cosh/C2 cosh If 3 



(67) 

(68a) 
(68b) 



a = 


-2zi|sin0|/|l-*ie w | 2 , 


(69a) 


b = 


(l-z 2 )/\l- Zl e M \ 2 , 


(69b) 


A = 


{a 2 + T] 2 + 2ar/costf)/|l - z x e^\ 2 , 


(69c) 


E = 


(2/? 7 sin l ?)/|l~z 1 e^| 2 , 


(69d) 


C = 


(2a?7sin^)/|l-z 1 e l ' 5 | 2 , 


(69e) 


F = 


(/3 2 + 7 2 + 2/37Cos^/|l-z ie 4,? | 2 , 


(69f) 


a = 


Zl + 2:2^3, 


(69g) 


P = 


z 2 + z 3 z\, 


(69h) 


7 = 


Z3 + Z\Z 2 , 


(69i) 


7] = 


-(1 + z 1 z 2 z 3 ). 


(69j) 



(Mlnll - Zl e v 



(70) 



Taking the limit L — > 00 and observing that A + > A_ , we 
obtain the bulk free energy per site (42) in the alternative 
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representation 



and therefore that 



f h (T)= lim f(°°\T,L) 

L— yoo 

= — ln(2 cosh Ki cosh K 2 cosh if 3 ) 



-/ M\n\l- z^] - — [ rfi91nA+. (71) 
!7r Jo 4tt J 



2tt 

Thus the excess free energy per site is 

f£°\T,L) = fg \T,L) + fg°\T,L), 

with 



(72) 



fT\T,L) 



N 



ln(cosh K 2 cosh K 3 ) 



1 r 2 ^ 

4ttWo 



and 



1 r 2 " 



di?ln 



p+ V A 



AT-1 



(73) 



(74) 



where we already anticipate that the latter expression 
does not contain any surface or interface terms. Near T c , 
dominant contributions to the integral in (74) arise only 
near $ = and $ = 2n. Since the integrand in (74) is 
invariant under 1? — > 2ir — it is sufficient to consider 
small positive 1?. Employing the small-t expansion (28), 
we obtain 



^(A-F) 2 +4\EC\ 

A-F-2a\E\ 
with 



A + F ^2(1- z 2 z 3 ) 2 , 

2(1 - z 2 z 3 ) 2 
N 

2(1 - z 2 z 3 ) 2 



\] X 2 + OJ 2 , 



N 



-x, 



(l-z 2 )(l-z 2 ) 
2(1 - z 2 z 3 )(z 2 + z 3 ) 



(75a) 
(75b) 

(75c) 
(76) 



and 



(l-z 2 )(l-z 2 ) Ndtf & Lt CooCl Lt 



l(l-z 2 z 3 )(z 2 +z 3 ) V m g) Z>U#l 



r sin 



- 1 cos 2 6»)(rcos 2 6» + r- 1 sin 2 6')] 1/2 5 



«[( 

= [l + i(r-r- 1 ) 2 sin 2 (20)]- 1/2 i, 



(77) 



where (8), (10b), (15b), (17), (18), (19b), (21), (32), and 
(33) have been used. It follows that 



A± w (1 - z 2 z 3 ) 



1± 



•s/x 2 + w 2 
TV 



P± 



\J x 2 + oj 2 ± a; 



2 v / 5 2 ~ 



(78a) 
(78b) 



A+ 2(1 - z 2 z 3 ) 2 



V / S 2 T 



P- 



\/x 2 + LJ 2 — X 



P+ Vx 2 + to 2 + x ' 
\ \ N ~ 1 

_ I ^ p -2V2 2 +w 2 

A+, 



(79a) 
(79b) 
(79c) 



From the expression for p + / A + we conclude that the sec- 
ond term on the right hand side of (73) receives a nonsin- 
gular contribution A r_1 [ln(l — z 2 z 3 )— \ ln2] plus a contri- 
bution to the integral that is not restricted to small 1? or 
2tt — 1?, since In^ + x/^/x 2 + uj 2 ) w al/u; for large w. This 
points towards a singular contribution that logarithmi- 
cally violates scaling of f^°\ and we will not consider 
this quantity in what follows. 

Combining (74), (76), and (79), we obtain in the large- 
N limit close to T c for the finite-size term 



(80) 



with 



lV\x) 



= -— [ dculn ( 1 
2tt J \ 



- x e -2^x 2 +LO* 



We note that 



j(/)(0)=/(+)(0)/4 = -7r/48, 



• (81) 



(82) 



which results from comparing (82) with (50) and (51). 

Combining (10b), (15b), (16), (19b), (32), (33), and 
(80), we obtain, close to T c , 



C\T,L)^^^f\x). 



(83) 



For free be, the free energy contribution /f s is dominated 
by its leading singular part /f S;S . With (1), (21), (35), 
and (77), we obtain for the asymptotic singular finite- 
size part of the free energy per bulk correlation volume 



C(^) = \ f „ (84) 
\l(i\x)/{Ax 2 ) for T<T C , 

where Eq. (77) gives x as a function of x, r, and 9. 
Eq. (84) represents a central result of the current sec- 
tion. Its most remarkable feature is that its right hand 
side depends only on the single variable x. We conjec- 
ture that the right hand side of (84) represents a universal 
scaling function in the bulk universality class of the two- 
dimensional Ising model. In Sec. IIIB2, we will explain 
its independence of r and 9 in the more general context 
of rf-dimensional films. 



7W)(S)/a 



for T > T c , 
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Next wc derive standard free energy scaling functions 
and critical amplitudes. Combining (16) and (18a) gives, 
in the asymptotic critical domain, 

Vco/£o = r sin 2 6 + r- 1 cos 2 6, (85) 

so that, multiplying (84) by L 2 /V co while observing (1), 
(21), (77), and (85), gives asymptotically 

L 2 C\t,L)=gW(x,r,6) 

= (rsm 2 6 + r- 1 cos 2 6)l( f \x), (86) 

with x from (77). We conjecture that C/( 00 ) is a universal 
function of x, r, and 9 in the bulk universality class of 
the two-dimensional Ising model. Note that the scaling 
function C/( 00 ) depends on the anisotropy parameters r 
and 9 only through the geometric factors on the right 
hand sides of (77) and (86). 

With (82), the related Casimir amplitude is 

= -^-{rsin 2 6 + r- 1 cos 2 6). (87) 
48 

For the isotropic case, Eq. (86) reduces to 

G™(x)=I«H<c), (88) 

and thus we may relate 

(oo) (£, r, 8) = (r sin 2 9 + r' 1 cos 2 9) G^\x), (89) 

with x from (77). With x = 2x, our result for Q[^\x) is 
identical to xiZ'°^ (x) for ordinary be provided at the be- 
ginning of Sec. 12.1.2 in [9] without derivation. Q^\x) 

is shown in Fig. 3. A derivation of from published 
scaling functions for the Casimir force is provided in Ap- 
pendix B 2. 

As another cross check, we compare our results with 
the Casimir amplitude found in Ref. [23], where it is de- 
fined by 

A INW = A 2 / fs , s (0,L), (90) 

which, observing (10b), (15b), and (32), is related to 
A(r,0)=0(O,r,0) by 

Ai NW = (diAi)A(r,0). (91) 

Ref. [23] provides the result 

A< N ° w = -(7r/4%, (92) 

with 

V = = ^(r sin 2 9 + 7--W0), (93) 

where Cj = cosh(2.Kj) and Si = sinh(2i\T i ). The second 
equality in (93) is obtained by using (16), (17), (18b), 



(19b), (25), (32), and (33), as well as the criticality con- 
dition 

SiS 2 + S 2 S 3 + = 1, (94) 

which is equivalent to (26) and (43). Comparison of (87) 
and (91)-(93) shows that our result for A<°°)(r, 9) agrees 
with the corresponding result in Ref. [23]. However, while 
A(°°)(r, 9) is universal, this is not the case for Ainw as 
defined in (90), since Ainw depends on lattice details 
as manifested by the factor d\jh\ appearing in (91) and 
(93). 

III. FILM IN d DIMENSIONS 

The results of Sec. II may be analyzed in the more 
general context of a c?-dimensional film of thickness L, 
i.e., for a system in an L x oo d_1 geometry near a d- 
dimensional bulk critical point. Note that, within this 
section, we refer without further mentioning always to 
the asymptotic critical domain. 

A. Bulk correlation lengths 

Weakly anisotropic systems may be related to corre- 
sponding isotropic systems by an anisotropic scale trans- 
formation, see, e.g., [29]. Such a transformation may be 
realized as a shear transformation (see [16, 18, 19] and 
[17] for the use of such a transformation in the context 
of critical phenomena for field-theoretic models and for 
<p 4 lattice models, respectively) 

x' = Mx, (95) 

where x and x' represent the d Cartesian coordinates of 
position vectors in the original anisotropic system and 
the related isotropic system, respectively, and where the 
d x d matrix M may be decomposed according to 

M = SR (96) 

into a rotation, represented by an orthogonal matrix R, 
and a subsequent rescaling, represented by a real diagonal 
matrix S. 

It is then straightforward to adapt the above transfor- 
mation to the case treated here, where the long-distance 
correlations of the anisotropic bulk system are described 
by a correlation lengths ellipsoid represented by a tensor 
E that is diagonalized by R, so that 

RER- 1 =diag(^,C 2 2 ,...,^), (97) 

with the correlation lengths £i,...,£<j along the prin- 
cipal axes of the ellipsoid, and where S is the volume- 
conserving (det S = 1) rescaling matrix 

S = Vl' d diagfo" 1 , £2" 1 , • • • , C d 1 ) , (98) 
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with the correlation volume defined by 

V co = ^ 2 ---U = (detS) 1 / 2 . (99) 
Thus the correlation length in the isotropic system is 

£ = v% d . (loo) 

Consequently, the squared correlation lengths tensor in 
the primed system is 

H' = MHM T = V 2 J d l, (101) 

with the d x d unit matrix 1. Note that 

ST 1 = V ( T 2,d M T M. (102) 

In the ip 4 lattice model and in the field-theoretic con- 
texts of Refs. [17] and [16, 18], respectively, the shear 
transformation is written as A _1 / 2 U, where U corre- 
sponds to our rotation R and A -1 / 2 corresponds to S 
with the qualification that, in contrast to our definition, 
A is not necessarily volume-conserving. In Refs. [16- 
18], a matrix A = U _1 AU was defined, so that A = 
A/(det A) 1 / 6 * may be used to parameterize the long- 
distance anisotropy [18]. Since we choose S to be volume- 
conserving, our conventions imply 

A = A= (M T M)- 1 = (dctH)- 1 /^, (103) 

with 

detA = detA = l. (104) 

Note that A is defined here through the physical cor- 
relation lengths in the asymptotic critical domain. In 
contrast, explicit versions of A were obtained in [16-18] 
by requiring the related shear transformation to lead to 
a transformed Hamiltonian, whose expansion in small 
wave numbers k is isotropic through order k 2 . While 
for standard (p 4 field theory these definitions should co- 
incide due to an exact mapping between the anisotropic 
and isotropic bulk Hamiltonians [16], there is no reason 
to believe that this procedure generally leads also for lat- 
tice models to A as defined here. We will return to this 
issue in Sec. IV. 

Let n be a unit vector that is orthogonal to the film 
boundaries and define the vector 

L = Lh. (105) 

For two-dimensional systems such as the Ising model 
treated in Sec. II, it is convenient to choose coordinates, 
where 

fi T = (l,0), (106) 

so that, with r and 9 as defined in Sec. II B, we have 

/costf sin6>\ 

sin 9 cos 9 J v ' 



and 




(108) 



from which follows with (96) and (103) 

i-i _ (r cos 2 + r- 1 sin 2 6 \ (r-r^ 1 ) sin(2fl) \ 
"* ~ \\(r-r- v )sm{W) r sin 2 9 + r" 1 cos 2 9) ' 

(109) 

Therefore, the parameterization of the anisotropy of a 
d-dimensional system by A and the corresponding film 
orientation by n reduces, for a two-dimensional system, 
naturally to the parameterization by r and 9, as mani- 
fested, for the two-dimensional Ising model, by Eqs. (59)- 
(61), (65), (86), and (87). 

B. Free energy 

For a d-dimcnsional system with restricted geometry, 
the system shape and the be are generally transformed in 
a nontrivial way by the shear transformation. A major 
simplification arises for films with (anti-)pcriodic, fixed 
or free be, since (i) a film is transformed into another film 
and (ii) the be are invariant. This means that we can ex- 
press film free energy scaling functions of the anisotropic 
system in terms of the corresponding scaling functions 
of the isotropic system with the same geometry and be. 
Then the modifications for the anisotropic system can be 
represented by geometric factors, as we already explicitly 
observed in Sec. II for the two-dimensional Ising model 
and as will be detailed for d-dimensional film systems 
below. 

For notational simplicity, we formulate what follows 
for the singular part f s of the free energy density /, even 
though for the particular case of the two-dimensional 
Ising model there is no unique separation of regular 
part f s and singular part / ns and consequently T does 
not exist, in contrast to A and/or Q, depending on the 
be. If universality continues to hold for anisotropic d- 
dimensional film systems, we expect Eq. (3) to be re- 
placed by 

L d f s (t,L)= F(x,A,n), (110) 

where J 7 is a universal function of its arguments and the 
imposed be. Eqs. (6) are replaced analogously. 

Due to the different ways a film transforms for periodic 
and antiperiodic be on the one hand and for fixed and 
free be on the other hand, we treat these cases separately. 



1. Periodic and antiperiodic be 

Consider a film with periodic or antiperiodic be and a 
length L of (anti-)periodicity in the direction represented 
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FIG. 4: Illustration of the shear transformation from an 
anisotropic to an isotropic film for d — 2, i.e., for an infi- 
nite strip. Shown are the behavior of the strip boundaries, 
the correlation ellipse, and the vector L under a rotation R 
and a subsequent rescaling S along the horizontal and verti- 
cal coordinate axes. L is the thickness of the resulting strip, 
needed for free and fixed be. For periodic and antiperiodic 
be, the indicated dots at opposite sides of the original, ro- 
tated, and rescaled strips, connected by the vectors L, RL, 
and L' = SRL, respectively, are physically identical. For 
these be, translational invariance allows to define new strip 
boundaries (dashed) that are perpendicular to the direction 
of (anti-)periodicity, so that L' = |L'| is the thickness of the 
resulting strip. 



by the unit vector n. The shear transformation M trans- 
forms the film into another film, albeit with a different 
vector 



L' = ML 



(111) 



describing the direction and length L' = |L'| of 
(anti-)periodicity of the isotropic system. For 
(anti-)periodic be, the orientation of the film boundaries 
is not unique. We choose them such that L is orthogonal 
to the boundaries of the original film and L' is orthogonal 
to the boundaries of the film with isotropic bulk correla- 
tion lengths, see Fig. 4. Other choices lead to skewed be. 
Results for such be are related to the results presented 
here by elementary geometric considerations. 

Now consider singular parts of free energy densities 
that exhibit scaling. For the isotropic model, we assume 



L /iso,s(^; L ) — J~[ so (x ), 



(112) 



with the scaling variable x' = (L'/^+.o) 1 ^- With the 
critical behavior (9) of the correlation length £l in the 
direction L of (anti-)periodicity of the original system, 
the scaling variables x from (8) and x' from above are 
identical, since with (101) and (111) we have 



(i/a) 2 = L T H _1 L 



L' T H' 



(L'/O 2 - (H3) 



Therefore, we may write for the free energy per correla- 
tion volume 



f s (t,L) = V co f s (t,L) = V co f iso , s (t,L') 

( x- dv T iso (x) for T > T c , 

~ \Rf{-x)- dv F iso {x) for T < T c , 



(114) 



where we have used Eq. (100) and that our shear trans- 
formation is volume- conserving and thus separately con- 
serves the correlation volume and the free energy den- 
sity. Analogous equations hold for / CX;S , /sf>, and /f s>s , 
if they exhibit scaling. Thus we find for d-dimensional 
films with periodic or antiperiodic be that the leading 
singular parts of (scaling) free energy densities per cor- 
relation volume depend only on the ratio of the length of 
(anti-)periodicity and the correlation length in the cor- 
responding direction and not on any other details of the 
shape or orientation of the correlation ellipsoid. For the 
two-dimensional Ising model this is reflected by the ex- 
plicit results (55)-(57). 
Multiplying (114) by 



L d /V co = {L/L'Y 



x dv for T >T C 

R d A-i) dv for T <T C 



and observing 

U = L|Mn| = L(h T A^n) 1 / 2 , 



(115) 



(116) 



which follows from (103), (105), and (111), we obtain the 
standard free energy scaling function 

L d f s (t,L)=T(x,A,h) 

= (n T A- 1 n)- d / 2 J- iso (i). (117) 

That is, T depends only through a geometric fac- 
tor on the relative orientations of the direction of 
(anti-)periodicity and the correlation ellipsoid. 

Specializing to d = 2-dimcnsional systems with conven- 
tions as in (106)-(109), Eq. (117) reads 

L 2 .f s (t,L)= T{x,r,6) 

= (r cos 2 6 + r" 1 sin 2 0) -1 .Fi S o(£). (118) 

Analogous equations hold for / ex ,s, ffs,s, an d / s f, s , if they 
exhibit scaling. This explains, from a more general point 
of view, the anisotropy dependence of the Ising model 
scaling functions presented in Eqs. (59)-(64). 



2. Fixed and free be 

Here we repeat the considerations of the preceding sec- 
tion for a film with fixed and/or free be, such as 00, ++, 
0+, or +- be. The treatment immediately extends to 
similar other invariant be. For definiteness, we choose 
coordinates where the "1" direction is normal to the film 
boundaries, i.e., where 



n 



(1,0,. ..,0). 



(119) 



As opposed to the case of periodic or antiperiodic be, the 
thickness of the isotropic film is no longer given by the 
length of the transformed vector L', see Fig 4. Instead, 
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according to Eq. (2.48) of Ref. [30], the film thickness L 
of the isotropic system is given by 



L = {dctA- 1 /dct[[A- 1 ]]) 1 / 2 L 1 



(120) 



where [[A -1 ]] is the (d—1) x (d— 1) right lower part of 
A -1 and where, with the conventions employed here, 
det A -1 = 1, compare Eq. (104). 
For the isotropic model, we assume 



L /iso,s(^7 ^) — J~i so {x), 



(121) 



with the scaling variable x = (L/^o) 1 ^*- From (100), 
(103)-(105), (119), and (120), we obtain 

(L/a) 2 = I/S^L = (A- 1 )!! detllA- 1 ]]^) 2 . (122) 

With x from above and x from (8) this translates to the 
relation 

x = {(A' 1 )!! det^A- 1 ]]}- 1 /^. (123) 

For the singular part of the free energy per correlation 
volume, we may write 

f s (t,L) = V co f s (t,L) = V co f is0tS (t,L) 

f x- dv F iso {x) for T > T c 



R^ d (-xy d ^ iso (x) for T<T C 



where we have used Eq. (100) and that our shear trans- 
formation is volume-conserving and thus separately con- 
serves the correlation volume and the free energy den- 
sity. Analogous equations hold for / 0X , S 7 /sf,s, and /f s , s , 
if they exhibit scaling. Note that as for the case of 
(anti-)periodic be, the free energy per correlation volume 
depends only on one suitably chosen scaling variable. For 
the two-dimensional Ising model with free be, this is re- 
flected by the explicit result (84). 
Multiplying (124) by 

and observing (104) and (120), we obtain the standard 
free energy scaling function 

L d f s (t,L)=F(x,A,h) 

= (dctllA- 1 ]])^^), (126) 

with x from (123). Note that the right hand sides of 
Eqs. (120), (122), (123), and (126) change their form for 
a choice of coordinates that does not imply Eq. (119). 
The general form may be obtained by replacing 

(A- 1 )!! -> n T A- 1 n, (127a) 
det(A~ 1 - nn T A^ 1 - A- 1 nn T ) 



det[[A-']] 



r^A-ifi 



(127b) 



Specializing to d — 2-dimcnsional systems with conven- 
tions as in (106)-(109), Eq. (126) reads 

L 2 f s (t,L)=F(x,r,9) 



(rsin^ + r^cos^-Fisotx), (128) 



with 



= [l + i(r-r- 1 ) 2 sin 2 (20)] 



-V(2") - 
X, 



(129) 



as obtained by combining (A 1 )n and dct([[A - 1 ]]) — 
(A-% from' (109) with (123). With v from (21), Eq. 
(129) reduces to the relation (77) found for the two- 
dimensional Ising model. Equations analogous to (128) 
hold for / CX ,S7 /sf.s, and /f SjS , if they exhibit scaling. As 
for periodic and antiperiodic be, the d-dimensional point 
of view provides an explanation for the anisotropy de- 
pendence of the Ising model scaling function for free be 
presented in (89). 

In Appendix B2, we sketch the derivation of the 
isotropic Ising model results 



Gl + o +) (x)=GL 0) (-t)=I {f) (-^ (130a) 



' ISO 

G^\x)=G^{2x)/A = I^{2x)/A, 



(130b) 



with and 1^ ' from (81) and (50), respectively. With 
( 124 ) x = 2x, our result for G^ +) (x) is identical to X^' + \x) 



provided at the beginning of Sec. 12.1.2 in [9] without 



derivation. G^ + ^ (x) and {x) are shown in Fig. 3. 
Together with (77), (124), and (128), we immediately 
obtain predictions for both the finite-size scaling behavior 
of the finite-size part of the free energy per correlation 
volume and the scaling functions of the finite-size part 
of the free energy for the anisotropic case for both ++ 
and 0+ be. Combining (128) with (130), while observing 
(51) and (82), the related Casimir amplitudes are 

A i++ \r,9) = -A(°+)(r,6>)/2 = A( 00 )(r, 0), (131) 

with A(°°)(r,0) from (87). With the definition (90), the 
results provided at the end of Sec. 2 in Ref. [23] for the 
Casimir amplitudes are, for ++, 0+, and 00 be, 



(o-l 



a (++) _ _A yu ^> II- IS 



,(0+) 



(00) 
INW 



(132) 



with A|° ^ from (92). While the results for a|++ } and 

^inw were derived in Ref. [23] only for rectangular lat- 
tices, the combination of (91), (131), and (132) shows 
that they remain correct for the general triangular lat- 
tice. 



IV. SUMMARY AND DISCUSSION 

We have investigated finite-size scaling of an 
anisotropic two-dimensional ferromagnetic Ising model 
on an infinite strip of width L with periodic, antiperi- 
odic, and free be in the direction perpendicular to the 
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direction of infinite extent. The model is realized on a 
triangular lattice with general couplings and lattice con- 
stants, so that the anisotropy may be realized by varying 
the lattice constants di and/or the couplings Ji. This 
allows for different microscopic realizations of identical 
anisotropic bulk long-distance correlations near the bulk 
critical point and therefore for a limited test of univer- 
sality. 

We find that the asymptotic scaling behavior of the 
finite-size contribution to the free energy per correlation 
volume may be described by a function of only one suit- 
ably defined scaling variable, cf. Eqs. (55) and (84). For 
periodic and antiperiodic be this holds additionally for 
the surface and excess free energy densities per correla- 
tion volume, cf. Eqs. (56) and (57), while, for free be, log- 
arithmic violations of scaling prevent the scaling behavior 
of the surface contribution to the free energy density. 

For periodic, antiperiodic, and free be, we provide ex- 
act scaling functions , Q^ a \ and of the finite-size 
contribution to the free energy density in Eqs. (59) and 
(86) and scaling functions and A^ of the excess free 
energy density in Eq. (61). We find that these functions 
only depend on variables related to long-distance corre- 
lations. We have chosen these long-distance variables to 
be the scaling variable x representing the ratio of the 
width L of the strip and the bulk correlation length in 
the corresponding direction (for T > T c ; for T < T c its 
analytic continuation is employed) as well as the aspect 
ratio r and the orientation angle 9 of the bulk correla- 
tion lengths ellipse. Consequently, the critical Casimir 
amplitudes A< a \ and A^ 00 ' provided in Eqs. (65) 

and (87) depend only on r and 9. Since our results for Q, 
A, and A are independent of the microscopic realization 
of the long-distance physics, we conjecture that they are 
universal functions of their respective arguments. That 
is, we expect identical functions to result for other mem- 
bers of the bulk two-dimensional Ising universality class 
for the same geometry and corresponding be. 

To understand the behaviors of the Ising systems de- 
scribed above, we have investigated general anisotropic 
systems on a d-dimensional film with the same be as 
investigated for the two-dimensional Ising model. We 
assume that the free energies under consideration ex- 
hibit scaling. Due to the choice of geometry and be, 
these systems have the simplifying property that their 
geometry and be are invariant under a shear transfor- 
mation relating them to a corresponding isotropic sys- 
tem. We find that the free energies per correlation vol- 
ume of such systems depend on only one suitably chosen 
scaling variable, cf. Eqs. (114) and (124). The relations 
between the scaling functions of the isotropic and the 
anisotropic systems are provided in Eqs. (117) and (126) 
and reproduce for d = 2 the relations found for the two- 
dimensional Ising model. For ++ and 0+ be, where no 
results for the triangular lattice Ising model are available, 
we point out that these relations, together with the ex- 
plicitly known isotropic results (130), lead to predictions 
for the corresponding anisotropic scaling functions of the 



finite-size contribution to the free energy density of the 
two-dimensional Ising model. We give explicit results for 
the corresponding critical Casimir amplitudes in (131). 

We have only treated geometries and be that are invari- 
ant under the shear transformation that relates the sys- 
tem under investigation to a system with isotropic bulk 
correlation lengths. In other cases we cannot expect to 
express scaling functions for the anisotropic case in terms 
of an isotropic-case scaling function with the same geom- 
etry and be. Rather, the isotropic scaling functions will 
depend on the geometry and be obtained by transforming 
the anisotropic to the isotropic case and the anisotropy 
will no longer be represented by mere geometric factors 
as in (59)-(61), (77), (86), (117), (123), and (126). How- 
ever, we still expect such scaling functions to be universal 
functions of their arguments. Within a given bulk univer- 
sality class, they should depend only on the scaling vari- 
able (s) used for the isotropic case, the asymptotic long- 
distance anisotropy represented by the matrix A from 
(103), the geometry, and the be. Hence we suggest that 
a quantity is universal, if it depends only on the be and 
on macroscopic physical observables such as the geometry 
and the macroscopic near-critical correlation lengths, but 
not on the particular microscopic realization from which 
the anisotropic long-distance critical behavior originates. 
Such a quantity should therefore be identical among the 
members of the bulk universality class under considera- 
tion. For practical measurements of, e.g., scaling func- 
tions, it is no longer sufficient to measure one correlation 
length as in the isotropic case. Rather, the measurement 
of correlation lengths in a sufficient number of directions 
is necessary to allow for the determination of the matrix 
A through (103). 

Let us put our arguments in perspective with the in- 
terpretations of [16-19]. From our arguments follows 
that it is not necessary to define universality only after 
a shear transformation to an isotropic system, as Diehl 
and Chamati suggest [19]. As noted in Sec. Ill A, Chen 
and Dohm [16-18] define their matrix A through micro- 
scopic parameters in the Hamiltonian so that, at least 
for lattice models, their matrix is only an approxima- 
tion to our definition which relates A through (103) to 
the asymptotic physical correlation lengths. Thus the 
dependence of scaling functions on their matrix A will, 
in general, not be universal, since it depends on the mi- 
croscopic realization of the anisotropy. In this sense it 
is correct when Chen and Dohm note that any depen- 
dence of physical quantities on their anisotropy matrix 
A is nonuniversal [16-18]. In contrast, we suggest here 
that quantities that are universal for the isotropic case, 
merely acquire an additional universal dependence on A 
from (103) for the anisotropic case. The parameters de- 
scribing A are called nonuniversal in [18]. This is correct 
in the sense that the relation of A from [16-18] to our 
A depends on the model under consideration. However, 
with our definition, A is merely an argument of, e.g., 
a scaling function, and should be viewed on the same 
level as the scaling variable x. It is then only the seal- 
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ing functions, whose universality can be tested and not 
that of their arguments x and A. Similar arguments 
hold for the parameters describing the shear transfor- 
mation, which are called nonuniversal in [18, 19]. This 
transformation merely describes a relation between dif- 
ferent physical situations (or different interpretations of 
the same statistical model), not necessarily one of them 
being isotropic. With our definitions, its classification as 
universal or nonuniversal is not a meaningful question. 

Our last argument concerns the validity of two-scale 
factor universality. Since A is a scale-free quantity, our 
universality interpretation does not interfere with two- 
scale factor universality in a formulation including the 
field h conjugate to the order parameter. For definite- 
ness, consider again film geometry with given be and as- 
sume that the free energy density exhibits scaling. Then 
we expect the universal isotropic scaling function [3, 5, 6] 



ergy in Sec. II C 1. 

On the one hand, we need, for sufficiently well-behaved 
functions /, the large-n results 



L d f s (t, h, L) = TidtL 1 '", C 2 hL^I»), 



(133) 



with nonuniversal scale factors C\ and Ci and critical 
exponents /3 and <5 to be generalized to an anisotropic 
universal scaling function 

L d f s (t, h, L) = F{C 1 tL y / v , C 2 hL f}S '\ A, n), (134) 

(L)—l/ v 

where the choice C\ = o makes the first argu- 
ment identical to x in (8) and leads to the h=0 limit 
J"(a;, 0, A, n) = F(x, A, n), compare (110). Since no new 
nonuniversal scale factor needed to be introduced into 
Eq. (134) as compared to the isotropic case, two-scale 
factor universality remains valid. We have not quanti- 
tatively considered this case, since no exact results for 
the free energy of the two-dimensional Ising model with 
a nonzero magnetic field are available. 

As an outlook, it would be interesting to investigate 
systems that do not exhibit the simplifying feature of 
both invariant geometry and be under shear transforma- 
tions. A possible example is the critical Binder cumulant 
for a two-dimensional ferromagnetic Ising model. Preci- 
sion results for this quantity are available for a square 
lattice in a square geometry with periodic be and a 45° 
anisotropy that is caused by a ferromagnetic [31, 32] or 
antifcrromagnetic [33] coupling on one of the lattice di- 
agonals, i.e., by different microscopic realizations. For 
recent results of such an investigation, see Ref. [34]. 
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Appendix A: Integrals and sums 

Here we collect some mathematical results needed for 
the determination of the scaling functions of the free en- 



f dxf(x)= J2 f(x,)S+^[f'(b)-f(a)]5" + 0(S 3 ), 

j=l/2 

(Ala) 

/ dxf(x) = £ f( Xj )S + -[/(&) - f(a)}6 

-^lf'(b)-f(a)}S 2 + 0(S 3 ), (Alb) 

respectively, where S = (b — a)/n and xj = a + jS. 

On the other hand, using the method of residues, we 
obtain in the large-n limit for real x 

- duj yjx 2 + Auj 2 - V ^Jx 2 + [(2fc + l)7r] 2 

= ^ +/(+)(*), (A2a) 

- / dco yjx 2 + 4u; 2 - V y/x 2 + (2/br) 2 

^-\+I { -Hxl (A2b) 

with J(±) defined in (50). 

The results (Al) and (A2) allow us to obtain Eq. (49) 
from Eq. (48). 

Appendix B: Isotropic scaling functions 

Here we relate the scaling functions Q and A to pub- 
lished scaling functions of the Casimir force for the 
isotropic two-dimensional Ising model in infinite-strip ge- 
ometry. This allows us to verify the isotropic limits of our 
results for and and to derive the scaling 

functions and G-^^ that can subsequently be used 
for predictions about the anisotropic case as explained in 
Sec. IIIB2. 

For the isotropic case, the Casimir force of a d- 
dimensional film system of thickness L is defined by 

d[L/ iso , ex (T,L)] 

-TCas.isoi-t , -k) = 



dL 

d[Lf iso4s (T,L)} 



(Bl) 



where the second equality reflects the fact that surface 
contributions do not contribute to the Casimir force. If 
its singular part -Fc as ,iso,s exhibits scaling, a correspond- 
ing scaling function A iso may be defined in the asymp- 
totic critical domain, so that 

L d F CaSiiso ^ s (t, L) = X iso (x) 

= (d-l)gUi)-v- 1 i dG} ^ 1 , (B2) 
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where (6b) and (8) have been used. With (21), this may 
be expressed, for the two-dimensional Ising model, as 



L 2 F Cas ,i so ,s(t,L) = Xi so (x) 



~2 d [Gi: 

-x 



,(£)/a 



dx 



1. Periodic and antiper iodic be 



(B3) 



for the isotropic two-dimensional Ising model in infinite- 
strip geometry to derive the scaling functions Giso\ 
and £; ( s o +) . From Eqs. (2.9) and (2.13) in Rcf. [39] 
follows for the isotropic-case Casimir force scaling func- 
tions for ++ and 00 be in our notation (X there is iden- 
tical to our x/2) 



Here we use results from the literature for X^J a ^ to 
verify the isotropic limits (62) of the scaling functions 
g(.p/ a ) provided in (59). Rescaling w — > \x\uj in (50), 
applying (B3) to the resulting expression for Q^ a \x) 
and subsequently scaling back according to uj — > to/\x\ 
gives 



du:\/ x 2 + 1 



[tanh(V^ 2 +w 2 /2) - 1] 
[coth(Vi 2 +w 2 /2) - 1] 



(B4) 



where the upper and lower expressions hold for periodic 
and antiperiodic be, respectively. Corresponding critical 
values are 



A-,r (o) - \ -;/ 12 ' 

I it/6. 



(B5) 



The upper result in (B4) agrees numerically with the 
Ising curve plotted in Fig. 15 of Ref. [35] (x there is iden- 
tical to our x). The upper and lower results in (B4) agree 
with Eqs. (3.16) and (3.19), respectively, of Ref. [36] (x 
there is identical with our x/2). The upper result in (B5) 
agrees with Eq. (1) in Rcf. [37], with Eq. (3) in Ref. [38], 
and with Eq. (58) in Ref. [27]. 



X^ + \x)=X\ 



(00) 
iso 



(-£) 



dw^/x 2 + w 2 



^ _|_ Vx 2 +ui 2 -x c 2j£ 2 +ui 2 
^/x 2 +tu 2 +x 



(B6) 



Equivalent results are obtained from Eq. (3) in Ref. [40] 
for both 00 and ++ be (x there is identical to our x) 
[41] and from Eq. (3.26) in Ref. [36] for 00 be (x there is 
identical to our x/2), while the results from Sec. 12.1.2 

of [9] for X^ a °\x) and X(^~\x) (x there is our 2x) are 
missing a factor 1/4. Using (B3) and observing that, 
for large positive or negative x, the scaling function Q 
contains by definition no surface or interface terms, i.e., 
no terms linear in x, we obtain from (B6) by elementary 

integration the scaling functions Gi^ provided in (88) 

and (130a) and G-^^ provided in (130a) . 

From Eq. (3) in Rcf. [40] (x there is identical to our 
x), we obtain [41] 



1 /"°° / \ 

^i ( s°o +) (») = ^ J dwy/x 2 +L0 2 (coth yjx 2 + lu 2 - lj 



4 a o ) (25)/4, 



(B7) 



2. Fixed and free be 

Here we use results from Refs. [39] and [40] for the 
Casimir force scaling functions X^\ X^ + \ and X^\ 



with Xfaf from (B4). With (B3) and taking again into 
account the absence of large- 1 x | linear terms in G, elemen- 
tary integration leads immediately to as provided 
in (130b). 
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